Spectroscopic Signatures of Massless Gap Opening in Graphene 
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Gap opening in graphene is usually discussed in terms of a semiconducting-like spectrum, where 
the appearance of a finite gap at the Dirac point is accompanied by a finite mass for the fermions. In 
this letter we propose a gap scenario from graphene which preserves the massless characters of the 
carriers. This approach explains recent spectroscopic measurements carried out in epitaxially-grown 
graphene, ranging from photoemission to optical transmission. 
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In the last few years the realization of free-standing sin- 
gle layers of graphene opened a new route to the possibil- 
ity of investigating the behavior of massless Dirac fcrmion 
in two dimensions. Indeed, in graphene the valence and 
conduction bands are formed by the p z orbitals of the car- 
bon atoms arranged on the two sublattices A and B of the 
honeycomb lattice. When the two sublattices are electro- 
statically equivalent the two bands meet at the K (K') 
points of the Brillouin zone, leading to a zero-gap semi- 
conductor with two conical bands = ±vpk (k — k|), 
which resemble relativistic Dirac carriers with zero mass 
and Fermi velocity vp . For this reason, K and K' are usu- 
ally referred to as Dirac points. Even though the massless 
Dirac spectrum makes graphene the perfect playground 
for investigating relativistic effects in quantum systems, 
for device application a tunable-gap semiconducting be- 
havior would be more suitable. Along this perspective, of 
remarkable interest are some recent experiments of An- 
gle Resolved PhotoEmission Spectroscopy (ARPES) in 
epitaxially-grown graphene, which reveal a finite band 
splitting 2 A ~ 0.26 eV at the K point P, 0. In these 
works the authors propose that the splitting arises from 
the inequivalence of the A and B sublattices, which in 
turn leads to a massive (ms) gapped spectrum [H, H| : 



±VM0 +a 2 



(i) 



The magnitude of the band-splitting gap has been re- 
ported to decrease by increasing the number of layers, 
or, from another perspective, by reducing the induced 
charge density. Similar ARPES spectra were reported 
previously in Ref. la , although a different interpreta- 
tion was proposed [a & 0]- The most convincing argu- 
ment in such a controversy comes from the electronic 
dispersion far from the Dirac point, which is at odd with 
the massive gapped spectrum of Eq. ([1]). Indeed, Eq. 
(HJ) predicts E™ s ± to be unaffected by the gap opening 
in the 3> A regime, and in particular the lin- 

ear asymptotic behavior of the upper band E™*\ m v-pk 
(^™+ ^ A) should match the linear behavior of the 
lower band Ej°Z. ~ —v-pk. A careful analysis of the 
ARPES data reveals on the contrary a finite off-shift of 



the two asymptotic linear behaviors [H, 0, [H, which 
cannot be explained even by the periodic modulation in- 
duced by the substrate 0] . Further discrepancies appear 
also in the profile of the dispersion at finite k x away from 
the K point (Ref. 1], Fig. S3): instead of the parabolic 
shape predicted by Eq. |T|), a "V" shaped conic profile 
appears, characteristic of a massless dispersion. 

Quite remarkably, the issue of the gap opening has 
been raised also by very recent optical-absorption mea- 
surements in epitaxial graphene[9]. Indeed, the analysis 
of the optical spectra based on the massive gap model ([T]) 
has two major drawbacks: from one side, it would sug- 
gest that no gap (or a negligible one) is present in the sys- 
tem, in contrast with ARPES results, from the other side 
it fails in reproducing the data in the visible frequency 
range, where the "universal" conductivity value of e 2 /Ah 
is expected theoretically EH EH, and tested experi- 
mentally in multi-layer graphite samples [l| or few-layers 
suspended graphene 14]. In this Letter, we propose a dif- 
ferent gap scenario that reconciles the gapped nature of 
the spectrum with the massless character of the fermions 
in graphene. As we shall see, such a gap model not only 
accounts very well for the ARPES spectra, but reconciles 
also optical measurements performed from the visible to 
the far-infrared (IR) regime. 

Let us start by introducing the Hamiltonian for free 
electrons in the graphene honeycomb lattice, in terms of 
the usual spinor tp^ — (c^jC^g). Linearizing around 
the K Dirac point (we put h — 1), and using k = 
fc(cos0,sin0), we can write: 
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whose eigenvalues are the usual gapless Dirac cones 
= ±v-pk. If the A and B sublattices are electrostati- 
cally inequivalent, an additional term oc Aa z (ai = i tXt y tZ 
being the Pauli matrices) should be added to H^, and 
the gapped spectrum of Eq. |T]) is recovered. In sim- 
ilar way, in Ref. 0] it was shown that a gap opening, 
related to inter-valley scattering, is always associated 
with a mass onset. Finally, off-diagonal intra- valley pro- 
cesses of the kind Qf a x + Q\a v were claimed to lead to 
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FIG. 1: (Color online). Electronic dispersion (a) and DOS 
(b) for the massless _E ml (solid black line) and massive i? ms 
(dashed red line) gapped models, with A = 0.13 eV. The 
dotted line refers to the Dirac-like dispersion e in the absence 
of any gap. 



a mere displacement of the Dirac cone off the K (K') 
points jij. However, this statement is correct only as far 
as a constant value of Qf , Q\ for k — > is assumed, 
while in general a momentum dependence of the off- 
diagonal self-energy cannot be ruled out. For instance, 
the (unscreened) Coulomb interaction leads to the off- 
diagonal self-energy Xk oc k log(/c c /fc)[cos <pa x + singer,,], 
where k c is a momentum cut-off for the Dirac-like coni- 
cal behavior [l5j. In this case, a correction to the linear 
Dirac dispersion is achieved, but no gap opens because 
X vanishes as k — ► 0. 

In this Letter we want to explore a somehow interme- 
diate possibility, where the off-diagonal self-energy has 
the structure: 



Xk — A [cos 4>a x + sin <f>a y ]. 



(3) 



It is straightforward to check that, in this case, a gap 
is opened at the K (K 7 ) points without the onset of a 
massive term. Indeed, Eq. (J3J) leads to the massless (ml) 
gapped spectrum of the form: 



with a corresponding density of states (DOS): 
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where Vbz = 5.24 A 2 is the volume of the Brillouin zone. 

While the main aim of this Letter is to focus on the 
phenomenological outcomes of this spectrum, a brief dis- 
cussion on the microscopic origin of the self-energy ([3|) 
will be given later. In Fig. Q]we compare the electronic 
dispersion (panel a) and the DOS (panel b) of the mod- 
els ((4]) and |T]). Two main striking features need here to 
be stressed: i) the gap in Eq. (jU, as mentioned, is cre- 
ated without affecting the conical electronic dispersion. 
In this sense no massive term is induced; it) in contrast 
to the case of the massive gap model (p}, the lower and 
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FIG. 2: (color online), (a) Intensity map of the spectral func- 
tion A(k,u) of the massless model [Eq. Q] for A = 0.13 eV, 
up = 6 eV Aand T taken from the experiments (see text). 
Dashed white lines show the mismatch between the asymp- 
totic linear behavior of the upper and lower bands, (b) Plot 
of the gap edge at finite k x for the massless and massive gap 
model (A,vf as in panel (a)). Dark squares are experimental 
data taken from Ref. [l|]. 



upper cones of the spectrum ([T]) are misaligned by the 
quantity 2 A. Such anomalous features are reflected in 
the DOS [Eq. (j^J)] where the linear vanishing of N(lu) 
as |cj| — > A + in the model ((4j) points out the absence of 
a quadratic massive term, while the misalignment of the 
cones is reflected in a corresponding mismatch of the lin- 
ear dependence of the DOS, which does not extrapolate 
to zero at the Dirac point. 

To make a comparison with the ARPES results in 
epitaxially-grown graphene, we show in Fig. [2^ the cal- 
culated spectral intensity J(k, ui) = A(k, oj)f(uo) of the 
massless model E± % where jj, is the chemical potential, 



f(u>) is the Fermi function and A(k, w) is the spectral 
function A(k, u) = -Tr{Im[(w+^+ir)/-.ff°-Xk]- 1 } = 
r/(27r){[(w -fi -v F k- A) 2 + T 2 ]- 1 + [(cj — fi+ v F k + 
A) 2 +r 2 ]- 1 }. In Fig.EHwe used \i = 0.4 eV, 2A = 0.26 eV 
(A = 0.13 eV)(H and we assumed a quasi-particle scat- 
tering rate Y[uS) = T + a\u>\, with T = 0.165 eV and 
a = 0.11 fitted from the ARPES dataQ away from the 
Dirac point using vf — 6 eV A[l6j]. Similar values could 
be estimated from Ref. [||. Notice that the presence 
of such a large scattering time at the Dirac point partly 
spoils the gap feature as k — > 0, because the two peaks 
of A(k = 0, ui) at oj — ±A overlap and cannot be re- 
solved, leading to the controversial interpretation of sim- 
ilar ARPES spectra in Refs. @,@| and respectively. 
The same feature is observed in Fig. [5^, where the ex- 
act dispersion at the Dirac point is not detectable due to 
the large electronic damping. Nevertheless, one can still 
clearly resolve the net misalignment between the upper 
and lower Dirac cones, which is peculiar to the model J3J), 
and it is in perfect agreement with all the existing experi- 
mental data [lib, 5]. Moreover, ARPES data taken away 
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from the K point provide us also with a direct evidence 
of the linear (massless) behavior of the electronic disper- 
sion at finite k. By measuring the dispersion at k y = 
and finite k x , as extracted from several k y cuts through 
the K point, one can easily discriminate between the two 
models ([T]) and (|3]). Indeed, within the model ^ the en- 
ergy spectrum a/A 2 + (vpk x ) 2 would appear parabolic 
within a momentum window k x < A/vp ~ 0.02 A -1 , 
while within the massless model (J4j) one expects a linear 
increase of the gap, A + Wi?|fc x |- The experimental data 
by Ref. [l| (Fig. S3 in Supplementary material) are re- 
ported in Fig. [2p: as one can see, the massive model E± s 
shows no resemblance with the data, while the massless 
model E± l allows one an excellent fit of the dispersion 
without any adjustable parameter. 

The massless character of the spectrum has also sig- 
nificant consequences on the structure of the optical con- 
ductivity <j(fl). To elucidate this issue we evaluate here 
the optical conductivity in the bare-bubble approxima- 
tion. As usual, a(fl) is given by two part s, associated 
to intraband and interband transitions [lfl [llj. At the 
leading order in r/|/z|,r/A <C 1, we obtain: 
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As one can see, the general structure of the optical con- 
ductivity is the same for all the gapped and ungapped 
models: a Drude peak of width r opt = 2r and an in- 
terband contribution which starts at a threshold given 
by the larger between 2fi and 2 A, and saturates at 

n 
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^> A, to the universal value e 2 /Ah [1 
When a gap opens, part of the spectral weight is trans- 
ferred to the interband transitions and the Drude peak 
decreases with respect to the ungapped case. In the 
model ©, <r(fi) can be obtained by replacing the fac- 
tors - A) — > (|/i| - A 2 /|^|) and (1 - 2A/|f>|)_-> 
(1 + 4A 2 /|fi| 2 ) respectively in Eqs. © and ([7]) (lb! Hi | . 
Two features allow thus one to differentiate the models 
(P) and ([3]): the relative weight of interband and intra- 
band contributions and the shape of the conductivity at 
the interband threshold. In particular, one can see that 
(i) the gap-induced reduction of the Drude peak is much 
stronger in the massless gap model (of order ~ A) than 
in the massive one (of order ~ A 2 /|/z|), see inset of Fig. 
GJi); (ii) while the massive gap model would give rise 
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FIG. 3: (Color online), (a) Optical conductivity a(Q) for 
(i = 0.1 eV, T = 300 K, V — 15 meV for the massless gap 
model (solid black line, A = 45 meV) and for the massive 
gap model (dashed red line, A = 73 meV). Inset: reduction 
of the zero- frequency optical conductivity in the two cases, as 
a function of A/fj,. (b): Comparison of the two models with 
the experimental data for T(lo) from Ref. 0. 



to a peak above the asymptotic value e 2 /Ah at the edge 
of the interband spectrum and to a rapid saturation to 
the asymptotic value, the massless model predicts a de- 
pletion in the correspondence of such edge (due to the 
factor 1 — 2 A/51 < 1), followed by a slow saturation to 
the universal value, see Fig. Such a depletion is 

again a consequence of the vanishing of the DOS at the 
gap in the massless model (see Fig. [lb), with consequent 
reduction of the spectral weight for transitions occurring 
between the two bands. 

The present results open thus a new perspective in 
the analysis of the optical properties of graphene grown 
epitaxially on SiC substrates, that can be deduced from 
measurements of the optical transmission T(u>), which 
in first approximation is related to the real part of 
the single-layer optical conductivity as T(w) = [1 + 

\l2\ , where nsic is the 
refractive index of SiC and N is the number of lay- 
ers. In Ref. indeed, transmission data for few-layers 
graphene samples in the frequency range from the far- 
IR to the visible were analyzed within the context of the 
massive gap model described in Eq. ([T]) , and the gap was 
concluded to be negligible within the experimental accu- 
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racy A < r opt ~ 10 meV Q. Such a fit reproduces the 
data in the far-IR to mid-IR range, but fails in the vis- 
ible range, where, according to this fit, the data would 
then indicated a conductivity larger than the universal 
value e 2 /Ah. The failure of the fit follow from the fact 
that the size of optical transmission T(S1) at SI w 3 eV 
and for — > 0, and hence the size of the optical con- 
ductivity in the corresponding range, are found to be 
of the same magnitude. According to the previous dis- 
cussion, this could be achieved by a transfer of spectral 
weight from the Drude peak to the interband conductiv- 
ity, as due to a gap opening. However, to reproduce this 
feature within the massive gap model one would need 
A « l/il-y/l - 7rr/2|^| = 87 meV, where ~ 0.1 eV, 
r = r opt /2 « 15 meV are extracted from the interband 
edge and from the width of the Drude peak . With such 
large value, A » T, the massive gap would be clearly de- 
tectable as a sharp peak at the interband edge, which is 
instead absent in the data. For this reason, the authors 
of Ref. Q extracted a vanishing gap from the fit based on 
the model (p}, which then fails in reproducing the data 
in the visible range. 

Such ambiguity can be naturally solved within the con- 
text of the massless gap model where: i) a relatively 
smaller gap value is needed to make the low and high 
frequency optical conductivity of the same magnitude; 
ii) the opening of a massless gap does not give rise to 
any peak at the threshold of interband transitions but to 
a depletion of the conductivity with respect to the univer- 
sal value, and then to a slow and smooth crossover toward 
the high-frequency regime, as observed in the data. This 
picture is in very good agreement with the actual exper- 
imental measurements. In Fig. we show the best fit 
to one set of the experimental data of Ref. Q by using 
the massless gap model. We take |/z| =0.1 eV, T = 300 
K and r = 15 meV from the experiments themselves and 
we estimate A = 45 meV, N = 18. For comparison, 
the fit with the massive gap model, constrained to repro- 
duce the experimental values of T(u>) in the low and high 
frequency limit, would give A = 73 meV and N = 18, 
and would result in a clear peak (shoulder) at the inter- 
band edge. We would like to stress that the even though 
the absence on an interband peak can also be accounted 
for by a vanishing gap, the similar magnitude of the low 
and high-frequency values of the optical transmission is 
a clear indication of a reduced Drude height, and hence 
of the presence of a gap. 

So far, we have analyzed in details the consequences 
on ARPES and optical spectra of the massless gap in- 
duced by a self-energy as described in Eq. ©. In the 
last part of this Letter we discuss the possible physi- 
cal origin of such off-diagonal self-energy. To this aim 
we expand the electron Green's function and the self- 
energy in terms of their Pauli components where, for 
the case of equivalent carbon sublattice, we can ne- 
glect the <Xz term. In the noninteracting case the off- 



diagonal Green's function, for instance oc a x , has the 
simple form: G x (k, iuj n ) — —kcos(p/[(ico n + fi) 2 — Vpk 2 ]. 
At the Hartree-Fock level, one can quite generally write 
S k = T J2w,n ^( k " k')G(k, iw„)e l "" 0+ , which can have 
off-diagonal components depending on the form of the 
potential V(k — k') (see for example the case of un- 
screened Coulomb interaction discussed in Ref. [TBI]). 
While a completely momentum-independent interaction 
would thus lead to a vanishing off-diagonal contribution, 
due to the angular average, a small anisotropy in the 
scattering angle can lead to a self energy of the form © . 
For example, one can use V(k — k') = V9{\4> — <j>'\ — <p c ), 
where 4> c = ir corresponding to isotropic scattering. In 
this case, one can easily check that Sk has exactly the 
form with A = V 2 k 2 sin 2 </> c /87r 2 . Further investi- 
gation on the possible source of anisotropic scattering 
could thus shed more light on the microscopic origin of 
the proposed self-energy ([3]). In this context, it is worth 
noting that a possible role of the doping is suggested by 
the observed vanishing of both the gap and the band mis- 
alignment as a function of the charge concentration^^]. 

In summary, we propose a gapped model for graphene 
which allows one to reconcile the massless Dirac char- 
acter of the carriers with the effects related to a gap 
opening at the Dirac point. We show that both ARPES 
and optical-conductivity measurements give clear indica- 
tions of such massless gap opening in epitaxially-grown 
graphene. Since ARPES measurements are available only 
for epitaxially-grown graphene, we cannot rule out the 
possibility that such a gap opening is restricted to these 
systems. Recent tunnelin g da ta on epitaxially-grown [l?] 
and suspended graphene [l8j are not conclusive on this 
respect: however, our predictions can be further tested 
experimentally and, if confirmed, they would pose strin- 
gent constraints on the interaction mechanisms at play 
in graphene. 
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